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ABSTRACT: A method is described to model the dynamics of tapered 
axial bar. of various crow .action, ba.ed on the well-known Craig/Bampton 
component mode .ynthe.l. technique. Thi. element i. formed in term, of 
the .tatic constraint mode, and interface restrained normal modes. T h' 8 i 8 
in contrast with the finite element, a. implemented »n NASTRAN where the 
Interface restrained normal mode, are neglected. The.e normal mode, are 
in term, of Bessel function.. Re.toration of a few of the.e mode. lead, to 
higher accuracy with fewer generalized coordinate.. The proposed model, 
are hierarchical so that all lower order element matrices are embedded in 
higher order element matrices. The advantage, of thi. formulation compare 
to standard NASTRAN truss element formulation are demonstrated through 
simple numerical examples. 


1. Introduction: Tapered bars and beams have high strength to weight ratios as well 
as architectural advantages. They are frequently employed to model structures in 1 - 
verse applications, such as ship masts, turbine blades, chimney structures or complex 
frame constructions. NASA (Langley) has tested a truss structure which is made from 
tapered members to be used in space applications [1]. The technical literature on tapered 
beams is indeed vast with a long history [2-10]. Tapered beam finite elements are ei- 
ther simple elements (e.g., Lindberg[2], Rouch/Kao[3]) having two degrees of freedom at 
each end or higher order elements (e.g., Thomas/Dokumaci[4], To[5]) havmg more than 
four degrees of freedom. Ovunck[6], Avakian/Beskos[7], Gupta[8], Banerjee/Wilhams[9] 
and Spyrakos/Chen[10] have uBed frequency dependent finite elments in their analysis o 
tapered bars. Banerjee/Williams[9] have developed exact dynamic stiffness matrices for 
Bernoulli- Euler beams. However the approach in References [6-10] involves the unknown 
frequencies of the overall structure. The general framework developed by Engels[ll and 
applied in References [12-14] allows for the derivation of hierarchical finite elemen s or 
any type of structural element. This approach does not require the prior knowledge of 
system frequencies, thus overcoming the need for an iterative procedure to compute the 
structural response. In the present paper, a dynamic finite element model for a certain 
class of tapered bars with loads acting only in the axial direction is developed. The ele- 
ment matrices are presented in parametric form and can be easily extended to formulate 
the finite elments for a wide variety of tapered bars covering most practical cases e 
convergence properties of this dynamic finite element, when compared to the regular finite 
element, are examined using numerical examples. 


2. Assumed Modes Method: The Lagrangian elastic displacement vector e(x,y,z,t) 
for a generic element can be written as the sum of two separate displacement vectors e/ 
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and e. 


e = cj + e 


(i) 


where e, i. a quasi-static displacement vector due to the interface displacements q, and is 
expressed as a Unear combination of static constraint modes <h, 


e/ = <f>i 9i 


( 2 ) 


The second part e represents the remainder of the total displacement vector e. It is that 
Lr'T.wlThi! measured relative to e, by an absolute observer Clearly, . vanishes 
L the „ coordinate, and therefore can be expressed as a linear combination of assumed 
modes <f> which are restrained at those qj coordinates 


e = <i> q 


( 3 ) 


The vector < represent, a set of generalized coordinates to be determined as part of the 
^'c^rdimd^^^should^bT^reBMd^tlmr^tlmugh^re^tr^ned 'normal often 

and (3) into Eq. (1) yields 

1 4/ 


! = [<t>I <f>] { 


def 


4>q 


( 4 ) 


so that . is written in term, of 

!^t C th“^p r renUUon of . in E,. (4) is complete in the sense that any degree of accuracy 
is theoretically possible as long as enough 4 modes are added. 

In the standard consistent mass matrix approach, the elastic displacement vector e 
over the element is represented as a linear combination of interpolation . « shape funcbona 
In fact these shape functions are identical to the static constraint modes and therefor 

model^l'nrte^of^ta^ne^lert^g't^e^dispUcemem^opcop r^am^ Umited numiier 

a— 

subdivk^onof basic element, is necessary, thereby simplifying the finite element grid and 
?31 thTmodd is hierarchical and therefore ha. all the advantage, associated with this 
mope'ty™ In addition, these finite element models are directly based on the assumed 
modes method which provides a sound theoretical basis. 
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In Reference [11] it is shown that for a linear elastic material, the element kinetic and 
potential energies T and V can be written as 


T = 




(5) 


V = K = j (B#) t C(B*) dv 

E 


( 6 ) 


in which the matrix C is the material stiffness matrix and B contains the appropriate 
partial derivatives in x,y and z. Futhermore, the matrices M and K represent the mass 
and stiffness matrices of a generic finite element E. In partitioned form, 


M = 


M n 

Min 

K = 

Kn 

0 ' 

(7) 

MJn 

Mnn J 


0 

Knn . 



where 


Mn = £ <j>J<t>i dm, Min = J <f>J<l>dm, Mnn = J <f> <j> dm 

r K E 


( 8 ) 


ana t r — — 

K n = J ( B<h) T C{B<h)dV, , K nn = J (B<f>) T C(B<f>)dV (9) 

E E 

Note that the K in partition is always zero, which means that no stiffness coupling exists 
between qj and q. 

At this point, a few remarks are in order. First, the matrices Mu and Ku repre- 
sent the standard finite element consistent mass and stiffness matrices for the element E. 
The consistent mass matrix approach represents in fact a static condensation or Guyan 
reduction whereby all noninterface degrees of freedom are eliminated. Secondly, if the in- 
terface restrained normal modes are used for the columns of <j>, then the present approach 
is identical to the Craig/Bampton component mode synthesis procedure as applied to a 
finite element. It should be emphasized that the element E is generic. This means that 
the proposed approach is valid, at least in theory, for any type of element. In the present 
paper, this general procedure will be applied to the special case of the tapered bars. 


3. Tapered Bars: The hierarchical stiffness and mas6 matrices of the tapered bar are 

obtained by solving the governing equations of motion for displacement. Figure (1) repre- 

sents linear tapered bar ab with a straight centroidal axis and the directions of the principal 

axes being the same for all crossections. The cross sectional area A(x) is given by 

A{x) = A a (l + c jY ( 3 * * * * * * 10 ) 

where c = d b /d a - 1 and Ai,d { (i = a, b) denote the cross sectional area and the depth 

respectively, c > -1 otherwise the beam tapers to zero between its ends and L is the 

length of the bar. 
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Although the formulation is valid for any n > 0, many 
practical cases of tapered bars arise when n is one or two, see 
Figure (2). If the geometrical properties of the element at both 
ends are given, the shape function for n can be derived as 

_ log[Ab/Aq) 

” log(db/d a ) 

For bars of closed box or I-section of constant width and vary- 
ing depth, n is not an integer and will vary slightly from E- 
q. (10) at all x other than the two ends. But the deviations 
are usually within one percent of the exact values. 



Tapered Element 


4. Static ran. tr.int Mode.: Consider the axial bar element a. lUuatrated m g 
ure (1). The bar is assumed to undergo vibration along .ts own axt. and as a ngri ^body 
only move along that same axis. The interface displacement, are defined as *(‘) - «(°- > 
and „(() = „(£, t), i.e., „ = |?i T and the displacement vector e » constdered to have 

only one component, i.e., u. Eq. (4) is therefore written as 


U = [ <f>l <f>2 \qi + ^ 9 


( 12 ) 




Figure 2. Sample Cross Sections 
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The static constraint modes, <f>i can be computed from the differential equation for 
the axial deformation u at x from end a of the tapered bar, 


iK)- 


where E is Young’s modulus. Integration of Eq. (13) gives, 

EA( X )% = c, 

Substituting for A(x ) from Eq. (10) and setting ( = 1 + cf gives 

du C\L 1 


Integrating Eq. (15) again, 


u = 


di EA a ci 

C X L 


where 


EA a C 


1 


/( 0 + ^2 


forn/1 

= In i for n = 1 


(13) 


(14) 


(15) 


(16) 


(17) 


The appropriate boundary conditions for the computation of static constraint modes are 
given by 

qi = 1, 92 = 0 (ig) 

9l = 0, 92 = 1 


The resulting static constraint modes are 


. no - /(i + <0 , no- no 
* l ~ no- /(i+c)’ H ~ no - hi + «) 


(19) 


Note that these constraint modes are in fact the shape functions used in the stiffness matrix 
of the tapered bar. 

At this point, enough information exists to compute M/j and K/j from Eqs. (8-9). 
Indeed, for the tapered bar, the kinetic and potential energies T and V are given as 


t = \ }<*■•) (£)’*• v - 5 }**•> (s) 1 * 


( 20 ) 


Substituting Eq. (12) into Eq. (20) gives 


Mn = J pA(z)^f<M*, Ku = J EA(x)-^- -j^ dx 


( 21 ) 
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(22) 


When n = 1, the explicit expressions for Kn and M/j are defined as 


Kn 


EA a c [1-1' 
L ln(l + c) [-I 1 . 


Mu = 


mu 

m 2 i 


mu 

m 2 2 


where 


-» “ [- 21n2(1 + c) - 21n(1 +c)+e ’ + 2cl 

---■ = iJ^ [(cI+2c+2),n(1+c) ^ 2 ' 2cl 

m 22 = r [2(1 + c) 2 (ln 2 (l + c) - ln(l + c)) + c 2 + 2c] 

4cln (1 + c) 


The counterpart expression for Kn when n ^ 1 is given as 


Kn = 


EA a c(n - l)(l + c) 


n — 1 


l -l 
-l l 


(»# i) 


(1 + c)”- 1 - 1 

The stifness matrix partitions Kji and Mu for n = 2 are evaluated as 
EA a { 1 + c) 

&II = r 


1 -l' 

pA a L 

2 

1 +c 

-1 1 

is 

II 

o> 

1 + c 

2(1 + c) 2 


(23) 


(24) 


(25) 


Because of the similarity of the governing equations between axial bars and torsional 
shafts Eq. (25) can be used as the stiffness and consistent mass matrices for tapered shafts 
by replacing the variables A a , E with J 0 , G respectively, where J a is the polar second 
moment of area and G is the shear modulus of elasticity. If it is decided that no extra q 
coordinates are to be retained in Eq. (4), then the procedure can be terminated at this 

stage. 

S. Interface Restrained Normal Modes: An entire class of hierarchical models can 
now be created solely on the basis of choosing the interface restrained assumed modes. In 
this paper, the set of interface restrained normal modes is used. The normal modes and 
their corresponding frequencies are obtained from solving the eigenvalue problem associ- 
ated with the partial differential equation, 


d_ 

dx 



(26) 


subjected to the clamped-clamped boundary conditions. 

Substituting for A(x) from Eq. (10) and letting { = 1 + cx/L gives, 

d 2 u n du pL 2 d 2 v. 

d? + tdt = e? W 


(27) 
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For harmonic vibration, 


u((,t) = U(()sinu>t 


(28) 


where t denotes time and u> is the circular natural frequency. Eq. (27) is modified using 
Eq. (28) as 

d?U ndU . 2 pi? tt n (29) 


de + i di 

The solution of Eq. (29) when c > 0 is 




( 7(0 = {AJ^aO + 


(30) 


where J and V are Bessel functions of the first and second kind and a — y \/£- 

For the case of n = 2, imposing the clamped-clamped boundary conditions in Eq. (30) 
yields the characterstic equation for the tapered bar, 

(31) 


The solution of Eq. (31) is 


sin ac = 0 


iir E . _ 

u> = —W — , 1 = 1,2,. ..,oo 

L\p 


(32) 


It is noteworthy that the natural frequencies of the tapered bar with clamped ends for the 
case of n = 2 are independent of c and are in fact the same as that of uniform rods. From 
Eq. (30), the interface restrained normal modes for n = 2 becomes 


- ^ sin («£ - a ) 
<Pi = C,- 

where the mass normalization constant C,' must satisfy 


(33) 


1 + C 

/ 


pAwUm = i 


(34) 


Substituting Eq. (33) into Eq. (34), C, is determined as 


C; 


pA a L 


(35) 


From Eq. (8), the mass matrix partitions are evaluated as 

L ) ~ 

M nn = f pA{xj$ T $dx = J, M/yv = [mji] , mji = J pA{x)<f>j<t>idx 

o 0 

% — 1 } 2 j > • • y OO, J 1)^ 


(36) 
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with 


ma . = \ggg (i + c )(_i)<-H (37) 

*7T Xir 

When c is set to zero, Eqs. (25) and (36) reduce to the case of uniform bars. Appli- 
cation of Eq. (33) into Eq. (20) gives the Btiffness matrix partitions, 

Kin = 0, Knn = Diag- [“'i • • •] 

and u>,- is given by Eq. (32). It can be seen that the tapered bar element matrices consist 
of very simple terms. 


0. Demonstration Examples: 


6.1 T&per»H Cantilever: This first example is concerned 
with the cantilevered bar clamped at the right end, see Fig- 
ure (3). The pertinent structural parameters are £ = 30 x 
10® psi, pg = 0.2839 lb/in 3 ,A a = 1 in 2 ,A b = 4 in 2 and L = 
72 in. 

The characteristic equation of this cantilever is given by 



a cos ac + sin ac = 0 (39) 

Figure 3. Cantilever 

Tables (1) and (2) show the number of converged frequencies 
to within a given percentage relative error when compared to 
the theoretical frequencies from Eq. (39) for different model 
orders n. 

Note that the bar is subdivided into the requisite number of elements to arrive at the 
model order n in the standard finite element method as implemented in CSA-NASTRAN 
whereas in the hierarchical finite element model, the number of bar elements is always one 
and the requisite interface restrained normal modes are added to arrive at n. 

Table 1. Standard Finite Element Method 
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Table 2. Hierarchical Finite Element Method 


e 

n 

< i % 

< 5 % 

< 10 % 

8 

16 

24 

32 

40 

6 

14 

23 

31 

39 

7 

15 

23 

31 

39 

7 

15 

23 

31 

39 


8.2 Planar Truss: Next, consider the planar single-bay truss 
as illustrated in Figure (4). The parameters are the same as 
in Section (6.1). The horizontal bar has a uniform cross sec- 
tional area of 4 in 2 . For planar truss elements, the transverse 
inertia must be taken into account. This time, however, no 
‘exact solution’ for the frequencies of the structure exists. A 
reference solution was obtained in two different ways: (1) by 
constructing a highly refined standard finite element model (2) 
by retaining a large number of normal modes in the hierarchical 
model. Both models were refined to the point where no signif- 
icant change in the frequencies occured and both the models Figure 4. Planar Truss 
produced the same results. 

Table (3) lists the frequencies from the hierarchical finite element model when two 
normal modes per bar are added (i.e., n = 16) as compared to the reference solution. In 
order to achieve comparable results from standard finite element method, each bar has to 
be subdivided into five beam elements (i.e., n = 66) and consistent mass matrix has to be 
generated within NASTRAN. 

Table 3. Frequency Comparisons For Planar Truss 



Mode 

Number 

Reference 
Frequency (Hz) 

Computed 
Frequency (Hz) 

Error 

% 

2 

2.606E2 

2.606E2 

-1.123E-2 

3 

2.999E2 

2.999E2 

-1.703E-2 

5 

6.958E2 

6.961E2 

-4.763E-2 

7 

1.467E3 

1.468E3 

-1.837E-2 

9 

1.757E3 

1.762E3 

-2.970E-1 

11 

2.141E3 

2.172E3 

-1.432 

13 

2.889E3 

2.898E3 

-3.317E-1 

14 

3.049E3 

3.081E3 

-1.045 

16 

3.382E3 

3.505E3 

-3.622 
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7. Conclusions For the case of tapered bars, it has been demonstated that the modal 
synthesis approach, where substructures are assembled to form the overall structure, can 
be used at the element level itself. This approach has superior convergence characteristics 
and the advantages of hierarchical formulation. 
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